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1.  INTRODUCTION 


An  ongoing  cfTon  is  being  conducted  at  the  US  Army  Armament  Research  and  Develop¬ 
ment  Cbmmand,  Ballistic  Research  Laboratory  (USA  ARRADOOM/BRL)  to  improve  the 
methodologies  that  are  used  to  calculate  the  vultterability  of  military  vehicles  to  ballistic  threats. 
The  Internal  Point  Burst  model'  whose  distinguishing  feature  is  its  ability  to  separate  the  debili¬ 
tating  cfTccis  of  spall  from  those  of  the  main  penetrator  constitutes  the  most  advanced  metho¬ 
dology  being  developed  at  the  BRL  and  elsewhere.  In  such  a  methodology,  the  vulnerability  of 
a  target  v.^hicle  is  calculated  as  the  expected  value  of  vehicle  incapacitation  due  to  the  direct 
impacts  of  primary  pcneirators  and  any  associated  metal  debris  fragments. 

In  a  simplified  version  of  a  more  detailed  stochastic  representation  of  the  point-burst  vul¬ 
nerability  process^^,  the  probability  equations  defining  the  expected  value  of  incapacitation  per 
projectile  (either  primary  penetrator  or  secondary  metal-debris  fragment)  for  a  vehicle  com¬ 
ponent  exposed  to  a  ballistic  threat  are  formulated  as  definite  integrals  whose  basic  form  is 
given  by 

X  = 

.wi.bi  lro,wo,bo,g)G(r,,w,,b,)</r,r/wi</b,</rodwo<^bo,  (lA) 

oo  oe  oo  oo  oe  oo 

J*JJ*5'(ro,Wo,bo)dro</worfbo=l.  (IB) 

oo  oo  oo 

In  this  cxprc.s.sion,  the  boldface  letiars  {ro,wo,bo)  and  (rj.Wj.bi)  are  the  values  associated  with 
the  random  variables  (Ro,Wo,Bo)  and  (Ri,W),Bi),  respectively,  where  R  quantifies  the  location 
and  W  quantifies  the  direction  of  motion  of  a  projectile.  Ctorresporxlingly,  the  set  of  variables  B 
is  used  to  quantify  (cbaracierke)  some,  but  not  all,  of  the  remaining  significant  features  of  a 
projectile.  The  subscript  0  is  used  to  identify  the  variables  associated  with  a  penetrator  at  its 
origin  and  the  subscript  1  is  used  to  identify  the  variables  associated  with  a  penetrator  at  impact 
with  a  critical  component  (Figure  1)  in  the  vehicle.  The  quantities  dr,  dw,  and  db  are  the 
infinitesimal  hypervolumes  associated  with  r,  w,  and  b,  respectively. 

The  quantity  S(ro,wo,bo)  is  a  continuous  furuition  used  to  give  the  probability  density  that 
the  rarxiom  variables  a.ssociated  with  some  projectile  at  its  source  will  have  the  values 
(ro,Wg,bQ).  The  function/(r|,wi,b|  I  ro,wo,bo,g)  is  the  unixirmalized  cortditional  probability  den¬ 
sity  that  a  projectile  created  as  (ro,wg,bo)  will  perforate  a  banier  characterized  by  the  set  of  ran¬ 
dom  variables  G  having  a  set  of  values  g  and  impact  the  critical  component  as  (ri,W],bi).  This 
latter  function  differs  from  the  conventional  definition  of  a  probabili^  density  function  (PDF) 
in  that  it  does  not  normalize  to  unity,  but  has  a  normalization  which  for  each  set  of  values 
(ro,wo-bo,B)  lies  on  the  range  from  0  to  1,  that  is 

®  ^  /i*//(ri,wi,b,  lro,wo,bo,g)dr]dwidbi  ^  1.  (1C) 

oe  oo  «e 

This  quantity,  herein  identified  as  the  perforation  PDF,  is  continuous  and  finite  in  b)  for  those 

'  J.R.  Rapp.  3nd  F  T.  Crown,  *An  Asve&smem  or  Existing  Point-Burst  Models  of  Armored  Vehicle  Vulnerability,*  Ballis¬ 
tic  Research  Laboratory  Memorandum  Repon  No.  02963,  Oaober  1979,  (UNCLASSIFIED).  (AD  #B043965L) 

^W.E  Beverly,  *A  Detailed  Stochastic  Ballbtic  VulnerabOny  Model  Tor  Armored  Mflitary  Vehides,*  Journal  of  Ballis- 
tia,  Vol.  Ill,  No.  3,  1979. 

^W.a  Beverly,  "A  Stochastic  Representation  of  the  Vulneiabiliiy  of  a  Critical  Component  in  a  MUhaiy  Vehide  to  Metal 
Fragments,'  Submitted  to  the  Journal  of  Ballistics  for  Publicatioii. 
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Figure  1.  The  Mathemtiicsl  Notuioa  in  t  Simple  SpiU  Dompie 


cases  where  a  binh  projectile  must  perforate  a  barrier  before  impacting  a  critical  componem  in 
the  vehicle,  and  is  the  product  of  Dirac  delta  fUnctionalS*  in  b|  for  those  cases  where  the  pro¬ 
jectile  can  travel  unimpeded  to  impact  the  critical  component.  In  either  case,  the  perforation 
PDF  is  Dirac  delta  fUncdonais  in  for  projectiles  impacting  the  critical  oomponeia  and  is 
identically  aero  over  its  total  range  for  projectiles  which  either  fail  to  petfoiaie  the  barrier  or 
impact  the  critical  component.  The  function  (?(ri,W|.bt)  is  any  well-behaved  function  whose 
values  on  the  impacted  surface  of  a  critical  componem  give  the  average  incapacitation  caused  by 
ail  penetrators  whose  "usetr  tandom  variables  have  the  values  (ri,wi,bi).  In  these  integrals, 
S(ro,wo,bo)  is  assumed  to  vanish  at  large  distances  frmn  its  maximum  value  so  that  an  integra¬ 
tion  over  all  values  of  all  variables  (identified  by  placing  the  infinity  symbol  at  the  bottom  of 
the  integral  sign)  will  yield  a  finite  expected  value  X  for  the  incapacitatioa 

•  The  dimensionality  of  these  integrals  and  the  complexity  of  the  geometrical  and  composi¬ 
tional  features  of  the  target  vehicles  from  which  these  integrals  are  deriv«d  could  make  their 
evaluation  by  deterministic  methods  prohibitively  expensive.  An  alternate  Monte  Carlo  method 
has  been  outlined  by  Reverly^  which  could  greatly  increase  the  efficiency  of  vulnerability  calcu¬ 
lations.  The  objective  of  this  study  is  to  artalyze  the  Mome  Carlo  procedures  used  in  Reference 
5  and  to  illustrate  their  use  by  evaluating  simple  definite  integrals 

In  the  next  section,  we  will  initially  analyze  the  Mome  Carlo  evaluation  of  simple  one¬ 
dimensional  integrals  We  will  then  extend  the  procedures  developed  for  the  one-dimensional 
case  to  multiple  integrals  We  will  also  analyze  integrals  having  the  form  illustrated  in  equation 
lA.  Then,  in  Section  ITT,  simple  imegrals  will  be  constructed  which  have  the  form  of  the 
integrals  discussed  in  Section  IT.  These  imegrals  will  be  solved  using  the  outlined  Monte  Carlo 
procedures  and  the  results  will  be  compared  with  the  known  closed  form  solutions. 


II.  THE  EVALUATION  OF  DEFINITE  INTEGRALS  BY  USING 
THE  MONTE  CARLO  METHOD 


A.  A  Simple  One-Dimensional  Intearal 

A  simple  one-dimensional  imegral  having  the  form 

X-y«^(x)G(jc)dx,  (2) 

where  //(x)  and  G(x)  are  well-behaved  (Unctions  on  the  interval  from  X)  to  Xj  ,  is  generally 
regarded  as  the  area  under  the  curve  H(jc)Gix)  from  Xj  to  xi.  Hovrever,  the  iruegral  can  be 
viewed  from  a  different  pengiecdve  if  the  integrand  is  rearnnged  to  obtain 

j^_^y//(x)G(x)  ax~cfH(X)f(jc)dx,  (3A) 

X,  L  _ 


where 


C 


— J^G(x)dx, 
*1 


(3B) 


*C.  Eisenhan.  and  M.  Zelen,  ’Elements  of  PfotaabOitr.'  Handbook  of  Physics.  EU.  Condon,  and  H.  Odishaw,  Ednots, 
McGraw-Hill  Book  Company,  Inc.,  New  York,  I9SB. 

^W.a  Beverly,  *The  Application  of  the  Monte  Carlo  Method  to  the  Solution  of  the  Internal  Pobtt  Bum  Vehicle  Ballis¬ 
tic  Vulnerability  Model,’  Ballistic  Iteseanch  Laboratory  Technical  Repon  in  Preparation. 
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and 

(3C) 

when 

X,  <X  <Xj, 

(3D) 

and 

/(x)-0. 

(3E) 

when 

X  <Xi,  X  >X2. 

(3H) 

The  integral  can  now  be  regarded  as  C  times  the  expected  value  of  Hix)  where  the  probsbility 

density  of  AVlues  for  x  is  giA«n  by  the  PDF/(x)  (Reference  4). 


Viewing  the  integral  from  the  second  perspecdMC  and  applying  the  strong  law  of  large 
numbers  (Reference  4),  the  integral  can  be  estimated  as 

(4A) 


where  the  Xy  are  a  series  of  mutually  independent  values  for  x  whose  common  PDF  is  /(x). 
This  mean  value  estimate  is  idemifl^  by  attaching  a  bar^to  X.  A  measure  of  the  confidence 
level  of  an  estimate  is  taken  to  be  its  standard  deviation  SX,  that  is^ 


8X- 


y-i _ 


/(y-1) 


1 

T 


(4B) 


According  to  the  central  limit  theorem  (Reference  4),  this  measure  of  confidence  can  be  inter¬ 
preted  as  predictii^  that  approximately  68  percent  of  a  large  number  of  similar  estimates  of  X 
will  fall  within  ±8X  of  X. 


A  step-by-step  outline  of  a  Monte  Oulo  evaluation  of  X  is  given  below  and  illustrated  in 
Figure  2. 

1.  Pick  a  value  Xy  by  sampling  the  PDF/(x).  A  variety  of  methods  have  been  developed 
for  conducting  such  sampling^’^  An  efficient  method  for  the  case  where  fix)  is  given  in  histo- 
gramic  form  is  to  pick  each  Xy  by  solving  the  integral  equation^ 

j/(x)dx-[W(0.1]y.  (5) 

*1 

The  quantities  Ry[0,l]y  are  a  set  of  independent  random  lumbers  where  each  laixtom  number 
is  picked  with  equal  probability  on  the  range  of  0  to  1.  The  reader  should  note  that  the 


Been,  Imrod-joion  to  the  Theory  of  Error.*  Addlion-WeUw  ?ubliihins  Compeny,  Inc.,  Reeding,  Men.,  1962. 

^W.  OJlier,  R.  Nagel,  R.  Goldstein,  ?S.  Miitebntn,  and  MJi.  Kalos,  *A  Geometric  Description  Technique  tor  Com¬ 
puter  Analysis  of  Both  the  Nudeer  and  Conwmionel  VubienbiUty  or  Aimoied  Military  Vehides,’  Methematicel  Appli¬ 
cations  Group,  Inc.  Report  No.  MAGI-670t,  August  1967, 

*Z.D.  Ceshwell,  CJ.  Everett,  and  O.W.  Rechard,  *A  Ptactical  Manual  on  the  Monte  Carlo  Method  Tor  Random  Walk 
Problems,’  Los  Alamos  Sctemi&e  Laboratory  Repon  No.  LA-2U0,  December  I9S7. 

^.A.  Schicidar,  'The  Mome  Carlo  Method,*  Pergamon  Praw,  Long  Island  City.  New  York.  1966. 
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Figure  2.  The  Moms  Cario  Eatimation  of  a  Simple  Oi»-Dimeiis»nal  Tnncral 


9 


sampling  as  deactibed  by  equation  S  generally  requires  that  the  PDF  be  normalized  to  unity.  We 
will  assume  that  this  method  is  used  in  the  following  outlines  of  Monte  Carlo  procedures. 

2.  Calculate  the  evem  score  Xj  as  C  Hixj). 

3.  Accumulate  the  score  in  a  bin  reserved  for  this  opeiatioa 

4.  Calculate  the  square  of  the  preceding  score. 

5.  Accumulate  the  squared  score  in  another  bin  reserved  for  this  operatioa 

6.  Reiterate  steps  1-S  for  a  total  of  J  events. 

7.  Calculate  X  using  equation  4A. 

8.  Calculate  the  standard  deviation  BX  using  equation  4R. 

9.  Determine  if  the  confidence  level  of  X  as  determined  ^  BX  is  adequate.  Conduct  more 
events  and  merge  their  results  wth  those  already  obtained  if  BX  is  too  large. 

Step  9  completes  an  outline  of  a  Monte  Carlo  evaluation  of  a  one-dimensional  integraL 
The  procedures  outlined  above  will  be  extended  in  the  following  Subsection  to  the  evaluation 
of  multiple  integrals. 


B.  Multiple  Tntettrals 


A  multiple  integral  of  dimensionality  /  which  has  the  form 

X-/*  ••  J**  “  flUxu  ...  .X,,...  ,x/)* 

>«,  A,  A, 

G(xi, .  ..  ,Xi,.  ..  ,X/)dxi  •••«*;•••  Ok/, 


(6A) 


can  be  evaluated  by  reiterating  the  procedures  outlined  in  Section  ITA  for  evaluatir^  one- 
dimensional  integrals.  The  imegratiotis  are  assumed  to  be  started  at  the  innermost  integ^  and 
conducted  toward  the  outermost  integral  where  I  is  the  ruruiing  index  over  the  variables.  The 
A,  and  B,  define  the  boundaries  of  the  integration  region  accordingly,  that  is 


Ai  ’’constant, 
A2=  Aiix^), 

m  m 


A,  =>4,(xi, . 


.  .  .Xi_i\ 


0 

0 


Bi  ‘’•constant'. 


•  • 

Bf^B^txi, .  . . 


•  • 


Ai  •’’AfCxj, . 


I  Xj *1 3,  B^  ^’Bf^x^ , ...  ,  Xf , ,  xj  _  j 


(6B) 


We  will  first  rearrange  the  integrarKl  in  equation  6A  to  obtain  : 
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/•••/•••  J/fOfi. ...  .x,,...  ,*,)• 

/(jfi . X, . x,)A\  ■  ■  ■  dx,  •  dXf, 


where 


Cl  -/*••/••  •  /g(x,.  . . .  .xi. .  .  .  .x,)dr,  ■  A;  •  •  •  dTf, 

/<)  ^ 

Gixu  ...  ,Xi.. ..  .X/) 


and 


/(xi. .  . .  .j;, .  . .  ,x/)- 


(7A) 

(7B) 


(70 


when  die  x;  have  values  located  within  the  r^un  bounded  by  the  and  B,  ,  and 

/(Xi,...,X, . X;)— 0,  (TD) 

when  the  x;  have  values  which  lie  out^de  the  region  bounded  by  the  A,  and  Bi.  Then,  using 
the  same  perspecdve  as  that  used  in  Subaacdon  ITA,  the  muldple  integral  of  equation  7B  can  be 
interpreted  as  the  expected  value  of  Hixt. . . .  ,Xt, . .  .  ,Xf)  where  the  variables  x;  ha>«  values 
predicted  by  the  joint  PDF/(xi,  .  .  .  .  . .  ,Xf). 


Howe\«r,  in  a  Monte  Carlo  esdmadon  of  X,  a  sample  value  (x;  )j  for  each  x,  during  trial  J 
has  to  be  piclced  individually  and  in  its  turn  The  probability  density  of  values  for  Xj  is  the  mar¬ 
ginal  PDF  /i(xi)  which  is  obtained  by  the  integration 
*1  ^ 

/i(xi)— J  J/(xi, . . .  ,Xi, . . .  .X/yOxi  •  •  dXi  •  •  -  dxt.  (8a) 

/<}  4,  A, 


Now,  similar  to  the  procedure  used  in  Secdon  ITA,  a  sample  value  (x\)j  can  be  derived  from 
the  integral  equadon 


Cotuinuing,  a  marginal  PDF/2C(xi)y,X2]  is  then  constructed  for  the  variable  X2  as 

/2[(Xi)y,X2]- 

*  • '  J  ••*//[(*!)>. *2. . . .  . x,]<*cj, . . .  .dx,.. . .  ,*tf, 


(80 


where 

B,  B,  Bf 

Cz-J  •••/■"  J/[(-*i);.-Xi. .  .  .  .X>. .  .  .  .XfldXi  •  •  ■  Xi  •  •  •  X,.  (8D) 

A -2  A,  A, 

Similarly,  a  sample  value  (X2)/  is  dern«d  for  the  variable  x;  from  the  integral  equation 

(xp 

S  f2[(.x,)j,xi]dXi~mO,\]v,  (8E) 

where  Aj  now  has  the  constam  value 

^2-A2[(Xi);1.  (8F) 

This  procedure  is  condnued  uiuil  a  sample  value  has  been  picked  for  each  variable  xi .  A  score 
Xy  is  then  calculated  for  the  event  as 

11 


. . .  ,(x/)>1. 


(9) 


The  remaining  procedure  is  identical  to  that  already  outlined  in  Subsection  ITA  for 
evaluating  one-dimensional  integrals.  This  procedure  is  outlined  below  for  mult^le  integrals 
and  illustrated  in  Figure  3. 

1.  Pick  a  set  of  sample  values  for  die  variables  X/.  These  operations  are  performed 
in  the  following  steps: 

A.  Consmict  the  marginal  PDF  /{(X])  by  using  equanon  8A. 

R  Pick  a  sample  value  (x))y  by  sampling /](X])  (equatbn  8R). 

C  Construct  the  marginal  PDF/^Cxj)  by  using  equations  8C  and  8D. 

D.  Pick  a  sample  value  (xi)^  by  samplir^ /iOtj)  (equation  8E). 

£  In  a  similar  manner,  pick  a  sample  value  (x})y  for  each  remaining  x;  by  using  the 
appropriate  marginal  PDF/;(x, ).  The  reader  should  note  that  the  marginal  PDF  in  each  case  is 
deriv^  from  a  joint  distribution  of  dimensionality  decremented  by  one  from  the  preceding 
sampling  evem. 

2.  Calculate  the  event  score  Xy  as  Ciiff[(xi)y, . .  .  .(x;  )y, .  .  .  ,(x/)yl. 

3-9.  Calculate  X  and  8X  as  described  earlier  In  this  section.  These  steps  are  identical  to 
those  already  described  in  Subsection  TTA  for  a  one-dimensional  integral 

Step  9  completes  the  outline  of  a  Monte  Carlo  evaluation  of  multiple  integrals.  The  pro¬ 
cedures  outlined  in  Subsections  TIA  and  ITB  snll  be  applied  in  the  following  subsection  to  evalu¬ 
ate  a  multiple  integral  having  the  form  of  the  vulnerability  integral  illustrated  in  equanon  1. 


C.  The  Vulnerability  Iiuegral 

The  vulnerability  integral  of  equation  lA  can  often  be  evaluated  by  more  expeditious  cal¬ 
culations  than  those  usually  needed  to  evaluate  multiple  integrals  of  equivalent  dimensionality. 
The  division  of  the  integr^  into  a  source  term,  a  perfbration  PDF,  and  a  critical  component 
incapacitation  fbnaion,  where  the  source  lenn  is  a  function  only  of  penetrator  birth  variables, 
will  simplify  the  picking  of  sample  states  for  penetratora.  This  gain  in  calculational  efficiency  is 
even  more  pronounced  for  the  general  vulnerability  eqtration  where  the  integrals  similar  to 
those  used  in  equation  1  are  stacked  to  represent  the  different  stages  in  a  penetrator  history 
(Reference  5). 

We  will  analyze  the  Monte  Carlo  evaluation  of  the  vulnerability  integral  by  refoimulaung 
equation  1  A.  In  the  new  form, 

X- 

,wt,b|  lro,wo,bo,g)Q(ri,wi,bi)i/r]  rfwj  rfbj  rfrorfwodbo, 

«e  «a 

“J'J*J*'^^*’o»^Oibo)f?(ro,wo,bo)rfro«/worfbo,  (lOA) 


where 

^(ro,wo,bo)-X///(ri,wi,b,  lro,wo,bo,i)Q(r,,W|,b|)rfri  rfw|</bi.  (lOB) 

m  m  m 
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Figure  3.  The  Monte  Carlo  Estimation  of  a  Multi-Dimensional  Integral 
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If  Z>(ro,wo,bo)  were  a  known  function  of  the  variables  Cro>wo,bo),  then  k  could  be  evaluated  by 
using  the  Monte  Carlo  method  to  pick  sample  values  ftom  5(ro,wo,bo)  as  outlined  in  Subsec¬ 
tion  TTB.  However,  in  vulnerability  calculations,  Z>(ro,ws,bo)  is  not  generally  ■>  tractable  aid  an 
alternate  method  of  solution  must  be  useiL 


In  a  Monte  Carlo  solution  still  based  on  the  strong  law  of  large  numben,  the  multiple 
integrals  in  equation  lOA  could  be  evaluated  by  picking  sample  birth  projectiles 
l(ro)7,(wo)>,(bo)yl  from  the  source  term  ^(ro.wo,bo)  by  using  the  procedures  outlined  in  Sub¬ 
section  rm,  and  then  evaluating  each  integral  i>[(ro>^,(wo)^,(bo)/]  which  is  associated  with  a 
sample  binh  projectile.  The  Mome  Carlo  evaluation  of  ea^  i>[(r^y,(wo)>,(bo)>l  could  also  be 
corduaed  by  using  the  procedures  outlined  in  Subsection  TTR  That  is,  a  sufficient  number  K 
of  sample  residual  projectiles  is  pidced  from  each  normalized  perforation  PDF 
/^fi,w,.bi  t(ro);,(wo)>,(ho)y.gyI,  defined  below,  to  obtain  an  accurate  estimate  Dj  of  the  asso¬ 
ciated  integral  £>l(ro)y,(wo)/,(io)/J  where  each  estimate  is  given  by 

Dj—{Pc)j  £0[(ri)/k,(wi)/k,(bi)jkl.  (IlA) 

k-l 

The  quantity  X  can  then  be  approximated  as  the  mean  of  a  large  number  of  Dj,  that  is 

Wj  ZC?[(rt)>.(w,)>k,(b,)jil.  (IIB) 


The  normalized  perforation  PDF  is  given  by 


/^[ri.wi.bi  l(ro)y,(wo)y,(bo)/,gyl 


/[ri.wi,bi  l(ro)y,(wo)^.(ho)^,g;l 

Wj 


(lie) 


where  the  rxirmalization  factor  (.Pc)y  is  the  probability  that  a  projsctile  birthed  as 
[(ro)>,(wo);,(bo)y)J  will  perforate  the  barrier  aid  impact  the  critical  component,  that  is 


Ff  I  (ro)y,(Wo)y,(bo)y  I  •“ 

,Wi,b]  I  (ro)y,(W())y,(bo)y,gy ]  r/Zi  dWi  db|.  (1  ID) 


Howev^,  the  approximation  implied  in  equation  1  IB  isn't  necessary  and  the  absolute  con- 
vergetce  of  X  to  X  can  be  obtained  ^  applying  the  Weak  Law  of  Large  Numbers  for  the  case 
where  TchebychefTs  Theorem  is  applicable  (Reference  4).  In  a  Monte  Carlo  estimation  based 
on  this  law,  X  is  estimated  as 

j^-7Z('’c)yO[(ri)y,(w,)y,(b,)yl  (12) 

■'y-i 

where  only  otK  sample  set  of  values  I(ri);,(wi)y,(b|  )y]is  picked  from  each  normalized  perfora¬ 
tion  PDF  /^[rj,w,,bi  ((rQ)y,(wo)y,(bo)J.  As  noted  in  Reference  4,  the  statement  concerning 
the  convergeiKe  of  the  sample  mean  X  to  the  universe  mean  X  is  weaker  when  only  the  weak 
law  of  large  numbers  is  ^piicable. 

A  step-by-step  outlitK  of  a  Monte  Carlo  estimate  of  the  vulnerability  integral  (equation 
1  a)  is  outlined  below  and  illustrated  in  Figure  4.  Stacked  vulnerability  integrals  such  as  those 
found  in  vulnerability  methodologies  (Reference  S)  can  be  estimated  by  reiterating  the  follow¬ 
ing  procedure. 

1.  Pick  a  birth  projecdle  by  sampling  5'(ro,wo,bo).  This  sample  projectile  is  identified  as 
[(ro)y,(wo)y,(bo)y].  Since  [ro,wo,bo]  are  each  assumed  to  be  composed  of  several  variables,  the 
procedures  outlio^  in  Subseaion  IIB  will  have  to  be  used. 
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START 


Figure  4.  The  Monte  Carlo  Estimation  of  the  Expected-Value-of-KiU  Vulnerability  Integral 
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1  Determine  the  values  gy  for  the  variables  associated  with  the  barrier  which  is  impacted 
by  the  projectile. 

3.  Derive  the  perforation  PDF  /[rt.wj.b]  l(ro)y,(wo)y,(b^)y,|y]  which  predicts  the 
description  (states)  of  the  residual  projectile  at  impact  with  the  critical  component 

4.  Determine  if  a  residual  penetrator  impacts  the  critical  component  by  comparing  a  ran¬ 
dom  number  RN(0,l)y  with  Set  the  event  score  to  zero  and  go  to  step  8 

if  a  residual  penetrator  does  not  impact  the  crid^  component  This  "Russian  Roulette*  played 
with  the  residual  fragmem  differs  from  the  procedures  implied  by  equation  12,  but  the  two 
methods  will  both  converge  toward  the  uue  value  of  X. 

5.  Pick  a  residual  penetrator  at  impact  with  the  critical  component  by  sampling  the  nor¬ 
malized  perforation  PDF  /^[ri.wi.bi  l(ro)y,(«o)>,(hD)y,gy].  This  sample  residual  penetrator  is 
idemifled  as  [(r,)y,(w,)y,(hi)y]. 


6.  Calculate  the  incapacitation  of  the  critical  componem  expeaed  from  the  impact  This 
incapacitation  score  is  identified  as  Xy  and  is  given  by 

X;-(?((r,)y.{w,)y.(b,)y].  (13B) 


7.  Accumulate  Xy  and  \}  in  bins  reserved  for  these  operations. 

8.  Conduct  a  total  of  J  similar  everus  by  reiterating  Steps  1-7. 

9.  Calculate  an  estimate  of  X  as 

X-if  Xy. 

•'y-1 


10.  Calculate  an  estimate  of  8X  as 


«X- 


;-i _ 


JU-l) 


(130 


(13D) 


11.  Assess  8X  to  determine  if  it  is  too  large.  If  necessary,  c^culate  more  histories  and 
merge  their  results  with  those  already  calculated  in  order  to  reduce  0X. 

Step  11  completes  the  outliie  of  a  Monte  Carlo  evaluation  of  a  multiple  integral  having 
the  form  illustrated  by  the  vulnerability  imegral  used  in  equation  1  A.  The  procedures  outlined 
above  will  be  applied  in  the  next  section  to  evaluate  a  definite  integral  constructed  by  using 
simple  analytic  functions. 


m.  SAMPLE  PROBLEMS 


A.  A  Monte  Carlo  Solution  Which  Uses  a  General  Multivariate  PDF 

We  have  devised  a  sample  two-dimensional  integral  which  can  be  used  to  illustrate  the 
procedures  outlined  in  both  Sections  ITB  and  ITC.  In  the  selected  definite  integral. 


-  225/16  - 14.0625, 


(14A) 


I  1 
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the  procedures  of  Subsection  TTB  are  illusinted  by  first  resnuging  the  integrand  to  the  form 

(14B) 

^11  » 

Equation  14B  can  then  be  further  chaagad  id  the  form  used  by  equation  7A,  that  is 


X—Ci  J*  J* Hbc^)/ixjf)dx4yt  (15a) 

by  taking: 


r 

0-4, 

(15B) 

/fCrj')— x^y*. 

(150 

and 

(15D) 

when 

(1,1)<(xo7)<(2,2), 

(15E) 

and 

/(xo»)-0, 

(15F) 

when 

(xo»)<(l,l),  (xj')>(2,2). 

(15G) 

The  constant  Cj  ■■9/4  is  inooduoed  to  omraalizB/Cxo')  id  uzdiy  over  the  integration  range  of 
the  integral,  that  is, 

/ //(jcor)* <fy dxdty-l.  (16) 


A  Monte  Carlo  estimate  of  A  can  now  be  caiculaad  by  picking  sample  'values  of  (xj^)  from 
/(xj')  and  then  calculating  the  aasodaied  loores  using  /TCxo')-  The  calculational  procedure  is 
given  below  and  is  illustrated  in  Hgute  S. 

1.  Pick  a  set  of  values  (.Xj^j)  by  sanqjling  /(x  j*)  owr  the  integration  range  of  the  pxob* 
lent  This  procedure  is  accomplished  in  the  following  steps; 

A.  Construct  the  marginal  PDF /((x)  u 

<17a) 

where  /\  (x )  is  taken  to  be  zero  when  x  lies  outside  the  integration  range  (1,2). 

R  Pick  a  sample  value  Xj  by  solving  the  integral  equation 

Xf 

J^<flC-^(0,l),y.  (17B) 

1  ^ 


C  Cbnstrua  a  marginal  POF/iCx^Ot)  fbr  as 

/jCXyOr)-— p — , 

V-J 


(170 
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Figure  5.  The  Mona  Carlo  Estimadon  of  the  Example  Problem 
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O.  Pick  a  sample  value  yj  for  y  bjr  solvios  the  integral  equation 


I 

The  quantity  (Cs)^  can  be  easily  shotm  to  be 

so  that  equation  17E  simplifles  to 

1  •> 


(17D) 


(17E) 


(17F) 


(17G) 


2.  Calculate  the  score  \j  for  the  sample  event  as 


(17H) 


1  Calculate  k/.  Accumulate  kj  and  k/  in  the  bins  reserved  for  this  operadoa 
4  Calculate  a  total  of  7  similar  simple  e\«nts  by  reitemting  steps  1'3. 

S.  An  estimaiB  X  of  X  is  calculated  as 

X— 7£Xy.  (ITT) 

•'  j-\ 


6.  An  estimatB  of  the  standud  deviation  SX  of  X  is  calculated  as 


ax- 


'Z.kj-Jk^ 

M _ 


7(7-1) 


(17J) 


7.  Assess  ax  to  deterraioB  if  it  is  too  large.  Tf  necessary,  calculate  more  hiSBuies  and 
merge  their  results  with  those  already  calculated  in  order  to  tedure  ax. 

Step  7  completes  an  outline  of  the  calculation  of  an  estimate  of  X  by  using  the  Monte 
Carlo  pfooedures  given  in  Section  TTR  A  compunr  progiam  MCTPl  for  perfbrming  this  calcu¬ 
lation  was  written  in  BASIC  and  is  giren  in  the  appendix  at  the  end  of  this  repoa  The  results 
of  cslcuiations  using  MCTP 1  ate  giren  in  TABLE  1  at  the  end  of  this  section 


R  A  Monte  Carlo  Solution  Which ' 


The  procedures  of  Subsection  TTC  are  iliustrated  by  first  rearranging  the  integrand  of  the 
example  inregial  (equation  14A)  to  the  form 


X  —  J* J*(x*)  (ay^)  (y )  ay  <dt 

I  I 


(ISA) 
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where  the  integration  is  assumed  to  proceed  from  the  inner  int^nl  to  the  outer  integral.  Equa¬ 
tion  18A  can  be  compared  with  equation  1 A  as 


X- 


S(sW(y\x)Q{y)tfydx~^jis^ix)/^(y\x)WOc)Qly)ifydx. 

i  1  1*1 

(18B) 

by  taking  ‘  j 

S(x)~-x^,  (19A)  ! 

e(y)-y.  (19B) 

1 

and 

i 

/(ylaf)=V-  (19C)  || 

The  quantity  W  (x),  identified  here  as  a  weighting  function,  is  introduced  to  compensate  for  1 

the  normalization  of  the  preceding  source  term  and  PDF  to  unity. 

We  define  a  normalized  source  term  5^(x}  as 

(19D) 

where 

1  '  1  1  ' 

(19E) 

Applying  the  Theorem  of  Bayes'®,  a  normalized  PDF/^(y  lx)  is  associated  with  by 
^_7x/t(ylx) 

1 

(19F)  j 

or 

(19G) 

Equation  18A  can  now  be  lewrinen  as 

1 

X  -  JJ  j  SHof)  Q(y)<<K 

(20A) 

and  rearranged  to 

1 

X=-J'JV(x)/f(y lx)  Q(y)4ydx. 

11  P  ^  J 

(20B)  1 

1 

The  weighting  funnion  is  given  by 


where  —  corresponds  to  the  quantity  /*<:  used  earlier  in  outlining  a  solution  of  the  vulnerabil¬ 
ity  integral  (Subsection  TTC). 

i 


Kim,  'StniMici]  Anilysh  fat  Mualon  and  Oacblon,*  ‘nia  Dndtn  Frast.  Ine,  Htradde,  OUnoh,  I97X 
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START 


Increase  J. 


<r  Isj-J  ?  >- 

_ 

Calculate  A,  and  SX . 
SA  too  large?''^ 

Print  Results. 

C  END  ") 


Figure  6.  The  Monte  Cario  Estimation  of  the  Example  Problem  where  the  Integrand  Is  Stnic- 
tured  in  a  Fbrm  Similar  to  That  of  the  Vulnerability  Integral 


Equation  20B  is  in  a  form  which  is  tractable  m  Monte  Carlo  evaluatioa  A  general  solution 
is  given  below  and  outlined  in  Figure  6.  The  steps  are: 

1.  Pick  a  sample  value  for  x  by  sampling  5^(x).  The  sample  value  is  deri'^cd  from  the 
integral  equation 


(21A) 


The  quantity  y  in  eqiuation  20A  is  sawd  and  used  lamr  as  a  factor  during  the  calculation  of 
the  evem  scoreL 

1  CbtKtrun  a  conditional  PDF/^(y  Ixy)  to  be  used  to  pick  a  sample  value  for  y.  The 
conditional  PDF  used  here  (equation  19G)  is  independent  of  x^,  but  this  independence  may  tM)t 
exist  for  cases  where  functions  other  than  are  used  in  the  construction. 

3.  Pick  a  sample  value  y^  by  sampling  f^iy  Ixy).  This  sample  value  is  derived  from  the 
iittegral  equation 


]>(y  lx,  )r<y 

1  f  * 


lx 

The  quantity  —  in  equation  20A  is  saved  and  used  later  as  a  factor  during  the  calculation  of 
the  event  score. 


4.  Calculate  the  score  for  the  evem  as 

5.  Calculate  X/.  Accumulate  K  and  X}  in  the  bins  reserved  for  this  operatioa 


(21C) 


6.  An  estimate  X  of  X  is  calculated  as 

_  t  ^ 

X-i^X,. 

•'  j-\ 


(21D) 


7.  An  estimate  8X  of  the  standard  deviation  of  X  is  calculated  as 

ixj-jl?  ^ 


1  Assess  8X  to  determine  if  it  is  too  large.  If  necessary,  conduct  more  histories  by 
reiterating  Steps  1-S  and  merge  these  results  with  those  obtairsd  earlier. 

Step  3  completes  the  outline  of  the  Monte  Carlo  evaluation  of  the  example  problem  by 
using  the  procedures  which  would  be  used  m  evaluate  a  vulnerability  imegr^  of  the  form 
shown  in  equation  1.  The  results  of  nst  calculations  of  this  example  problem  are  given  below 
in  Subsection  ITTC 
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C  Example  Probleni  ReauJu 


A  computer  program  was  written  in  BASIC  (APPENDIX)  to  solve  the  sample  inti^ral  hy 
using  both  Monte  Garlo  procedures  described  in  the  preceding  sections.  The  results  of  test  cal¬ 
culations  where  5000  histories  xrere  conducted  (br  each  estimate  are  given  below  in  Table  T. 
The  reader  should  note  that  standard  deviations  of  magnitude  ten  times  that  of  the  sample  cal¬ 
culations  could  be  obtained  by  using  approximately  fifty  histories.  Confidence  levels  of  this 
order  of  magnitude  would  usually  be  acceptable  in  vulnerabili^  calculations. 


IV.  CONCLUSION 


We  have  outlined  the  use  of  the  Monte  Carlo  method  for  solving  expected  value  integrals 
of  the  type  encountered  in  ballistic  vulnerability  calculations.  As  an  aid  to  the  comprehension  of 
these  methods,  the  outlined  procedures  (Subsection  ITB  and  nC)  were  applied  to  a  sample 
problem  for  which  a  closed  form  solution  existed  and  could  be  readily  calculated.  The 
insignificant  difTerences  between  the  computer-generated  Monte  Carlo  solutions  and  the  cx.ict 
solution  (Table  1)  serve  to  irxlicate  the  viability  of  the  Monte  Carlo  techniques. 

The  reasons  for  using  the  Monte  Carlo  method  in  preference  to  deterministic  methods  to 
solve  vulnerability  problems  can  be  summarized  as: 


1.  Regardless  of  the  dimensionality  of  die  integration,  a  Monte  Carlo  estimate  of  multiple 
integrals  converges  toward  the  true  value  as  where  N  is  the  number  of  trials  used  in  calcu¬ 


lating  the  estimate.  Therefore,  the  Monte  Carlo  method  may  be  the  most  efficient  procedure 
for  solving  expected-value  problems  when  the  dimensionality  of  the  associated  multiple 
integrals  is  large  and  confidence  levels  for  the  estimate  on  the  order  of  five  percent  are  accept¬ 
able. 


2.  The  probability  of  introducing  systemadc  errors  such  as  those  which  might  be  inuo- 
duced  in  deterministic  calculations  using  a  computerized  phantom  of  the  target  vehicle  is 
greatly  reduced. 

3.  Calculaiional  simplicity  can  often  be  obtained  by  the  division  of  the  imegrarxl  into  one 
part  used  in  calculating  scores  and  another  part  used  to  derive  the  sampling  PDFs.  In  parucu- 
lar,  the  PDFs  used  in  describing  the  ballistic  phenomena  are  often  directly  derivable  in  terms 
of  the  ’’used*  random  variables  B  and  G. 

4.  The  calculadons  can  be  organized  so  that  the  sample  events  correspond  in  a  one-to- 
one  manner  with  actual  ballistic  events.  This  organization  would  aid  vulnerability  analysts  in 
conducting  calculations  and  interpreting  results. 

5.  In  many  vulnerability  problems,  sample  events  can  be  obiaited  more  easily  than  the 
associated  PDF.  Such  problems  must  generally  be  solved  by  using  the  Monte  Carlo  method. 
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METHOD  OF 
SOLUTION 

RUN 

NO. 

ESTIMATE 

STANDARD 

DEVIATION 

VULNERABILITY 

INTEGRAL 

1 

14.058 

0.048 

2 

24.153 

0.046 

3 

14.043 

0.048 

MULTI-VARIATE 

SAMPLING 

1 

14.153 

0.102 

2 

14.062 

0.101 

3 

14.044 

0.100 

CLOSED  FORM 

14.0625 

Table  1.  The  Monte  Caiio  Estunates  of  Uie  &aniple-Piobleni  Integral 
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appendix  the  computer  prochiam  used  to  solve  the  example 

PROBLEMS 


i  PEM  Mt:  TEST  PROG  1  MCTPl 
iO  INPUT  "PI,.  N"  .  PI.  N 
;0  FOR  I=1T0  N 
ro  GOSUB  ■■  70 
40  X=SQP  l+i  +'R  '■ 

50  GOSUB  '70 
S  0  't*  “  Q  P  i  P  ■' 

70  Z=7!*V*X-+:V 
i.0  Si=Si+7 
'?0  S2=S2+Z'*-Z 

100  SELECT  PRINT  005:  PRINT  I,.  Z 

110  NEXT  I 

120  S1=S1.-'N 

120  S2=S2-N*S1*S1 

140  S  2 = S  2  /  '■'  N  ♦  N  - 1  ,■>  > 

150  S2=SOR''S2> 

1S0  SELECT  PRINT  215 
1'0  PRINT  5'*S1.'4,.  9>*<S2/4 
120  EfK' 

4000  C'EFFN  ■  70 

4001  PEN  PN  SP 
4005  P1=25+'P1 
4010  GOSUB  71 
4015  '■-1=5^P1 

4  0  20  G  0  SUB  71 
4025  P=P1.  S710£"5S4 
4020  RETURN 

4050  DEFFN ' 71 

4051  REM  RN  SP2 
4055  R2=Pl.  'S7lv3S8S4 
40K0  P”  =87108864*  I  NT  R2  'J 
4065  P1=R1-PZ 

4070  RETURN 


1  REM  MC  TEST  PROG  2  MCTP2 

10  INPUT  "Rl.  N"  .  Rl..  N 

20  FOR  I=lTO  N 

2-0  g0',^UE'  70 

40  X!=  *:  l4‘7*P  '0.  ZZZZZ3 

50  GOSUB  '70 

t0  V=  •:  1+7  +P  :■  '^0.  ZZZ33Z 

70  Z=X*V 

80  S1=S1+Z 

80  S2=S2+Z+'Z 

100  SELECT  PRINT  005:  PRINT  I 

110  NEXT  I 

120  S1=S1.  'N 

120  S2=S2-N*S1*S1 

140  S2=S2.'''  <  N*  <  N-1  >  > 

150  S2=SC!R''.S2> 

160  SELECT  PRINT  215 

170  PRINT  48*S1.'’9..  4S'*S2.-'S' 

180  END 

4000  DEFFN '70 

4001  REM  FT4  SP 
4005  R1=25+P1 
4010  GOSUB  -'Zl 
4015  F:l=5*Pl 
4020  GOSUB  71 
4025  R=F:l.-‘67108864 
4020  RETURN 

4050  DEFFN  71 

4051  REM  RN  SR2 
4055  R.2=R1. '67108864 
4060  R 2=67 108864*  I  NT  R2 
4065  R1=R1-R2 

4070  RETURN 


oisicuauo  pa(M  bumk-not  fiimu) 
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